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Abstract 

We give a precise definition for excitations consisting of a droplet of size 
n in the XXZ chain with various choices of boundary conditions, including 
kink boundary conditions and prove that, for each n, the droplet energies 
converge to a boundary condition independent value in the thermodynamic 
limit. We rigorously compute an explicit formula for this limiting value using 
the Bethe Ansatz. 
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1 Introduction 

In this paper we study the low-energy spectrum of the one-dimensional spin- 1/2 
ferromagnetic XXZ Heisenberg Hamiltonian in the thermodynamic limit. The spe- 
cific questions we are interested in concern the excitations that describe droplets, 
i.e., finite domains of reversed magnetization. The simplest case is where the in- 
finite chain is in one of its two translation-invariant ground states with all spins 
parallel or antiparallel to the z— axis. A droplet excitation is then a state with n, 
n > 1, opposite spins that form, up to quantum fluctuations, a compact cluster 
which moves through the system as a unit. Since the model also has kink and 
antikink ground states ^2 ^ in which two halfs of the chain have opposite mag- 
netization with a transition region in between, it is interesting to ask about droplet 
excitations with respect to such a ground state. This raises some interesting ques- 
tions about how to define a droplet excitation in this case and how to approximate 
them by excited states in finite volume obtained by imposing boundary conditions 
or a constraint such as a particular value of the total magnetization. 
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We will consider the spin- 1/2 XXZ chain of length L with Hamiltonian 

1 1 1 

x=l 

and study i?[i,L] +/ii,L, where hi^L is one of the following three choices of boundary 
term: 

periodic b.c: h^^L = [^1 Si + Sl Sl] - Sf Sl (1) 

droplet b.c: hi^L = S {Sf + Sl) (2) 

kink b.c: hi.L = -a {Sf - Sl) (3) 

See Section |21 for suitable choices of the constants a and 6 as well as other 
definitions. In the first two cases one can define the droplet energy by restricting 
the Hamiltonian to invariant subspaces of fixed total third component of the spin. 
Let A(n) denote the smallest eigenvalue of the Hamiltonian under consideration 
restricted to the subspace of states with n down spins and L — n up spins, < 
n < L, which is called the space of magnon states. For the case of periodic or 
droplet boundary conditions, the energy of a droplet of size n is then defined to be 
X{n) — A(0). For the model with kink boundary conditions this strategy does not 
work, since A(n) is attained in a kink ground state and is independent of n. The 
kink ground states form a multiplet of maximal SUq{2) spin: iSmax = L/2. It turns 
out that the correct subspace to define droplet excitations is the subspace of fixed 
total spin S = Smax—n. The mathematical explanation for this definition lies in the 
existence of a linear isomorphism between the space of n— magnon states and the 
"highest weight" vectors of "weight" S'max ~ n. The quotation marks are necessary 
here, since the isomorphism only exists for the infinite chain and the weights are 
not well-defined (infinite). See Section [^31 for the definition of this isomorphism, 
which we will denote by R. R is a bounded invertible operator that intertwines 
the Hamiltonians with kink and droplet boundary conditions on the infinite chain. 
Another, more physical interpretation, is that the total spin quantum number 
associated with a droplet of size n is n/2, in agreement with the case n = 1, more 
commonly known as spin waves [21 

Our main results can be summarized in words as follows: the droplet energies 
defined with the different boundary conditions above all converge to the same 
value in thermodynamic limit and that value can be computed exactly by the 
Bethe Ansatz. The result is given in Theorem 12 .11 In the proof of this theorem we 
use Perron-Frobenius type arguments to turn the Bethe Ansatz calculation into 
rigorous mathematics. For the precise definitions and mathematical statements we 
refer the reader to Section |21 

The main motivation for this study is to complete our understanding of the 
low-lying spectrum of the XXZ chain, which is important for a variety of problems 
involving the dynamics. As a by-product we have also come a step closer to a 
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complete proof of the completeness of the Bethe Ansatz in the thermodynamic 
limit. 



2 Set-up and main results 

2.1 The kink Hamiltonian 

For L G N+, consider a spin chain on the sites of [1, L] C Z. The Hilbert space is 
H. = = ®xe[iM '^x, where Hx is a two-dimensional Hilbert space for each 

X e We take an orthonormal basis of Hx to be the Ising basis {|t) , 

The spin-1/2 representation of SU(2) is defined on through the matrices 



1/2 
1/2 





i/2 



-i/2 




1/2 
-1/2 



in the {jf) , ||)} basis. For each x G [1, L] and i G {1, 2, 3} wc have the operators 
SI on Ti. where 5" acts on Hx and is tensored with 1 on Hy for ail y ^ x. 
The XXZ model is the Hamiltonian 



L-1 



H 



x=l 

1 



r»3 



[^x ^x+1 



^x ^x+l] 



(4) 
(5) 



with A > 1. Since we consider 1/A, it is allowable that A = +cx:. Let q be the 
number in [0, 1] such that A = {q+q~^)/2. Then a modification of this Hamiltonian 
is the so-called kink Hamiltonian 



H 



' c3 



a 



(6) 



where a is the constant 



a 



(7) 



1 + q^' 

When <7 — 1, this gives the isotropic Heisenberg model without boundary fields. It 
is useful to incorporate the alternating boundary fields into the nearest-neighbor 
interactions, 



1 



1 



''■x,x+l — 1^ ~ ^l^x+1 ~ IT {^l^l+l + '^x'^x+l) ~ 77 ('^x ~ '^x+l) • (8) 



4" ~x~x-t-i ^ 

Then the kink Hamiltonian can be written as 

L-1 
x=l 



(9) 
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When we want to emphasize the chain for the Hamiltonian, we will write i^fi^L] 
for H and Hf^ ^ for H^. 

There are three important operators commuting with each nearest-neighbor 
kink interaction, separately. The first is 

L 

x=l 

This is the usual total-magnetization operator for representations of SU(2). The 
other two operators are g- versions of the total raising and lowering operators 

5+.] :=i:'i-''''^-^'-''s^^ (11) 

x=l 

^[T,.] :-i:S-q'^''-^^-^'iK (12) 

x=l 

Note that these operators are only well-defined when < g < 1. These three 
operators together give a representation of the quantum group SU^(2). (For readers 
unfamiliar with quantum groups, we will present all the details necessary for our 
results.) 

The total magnetization eigenvalues are {L/2 — n : n = 0,1, . . . , L}. The 
eigenspace for the eigenvalue L/2 — n will be denoted as Ti,{n). It is an invariant 
subspace for H^. For n < L/2, S'^ ^ maps Ti,{n) isomorphically onto its image 
in T-L{n + 1). (For a proof of this and other facts about the representations of 
SUq(2) c.f. 0.) Moreover, this image is an invariant subspace of H^. Therefore, 
so is its orthogonal complement. For 1 < n < [L/2J, we may define 71^™ {n) as the 
subspace of 7Y(n) such that 

n{n) = 7^h*(n) © 5'[7.i]7^(ri - 1) . 

Define H''^(0) = H(0). Then, ^^^(n) consists of vectors in V. which have total 
SUq(2) spin equal to L/2 — n, and which are highest-weight vectors in the sense 
that S'jij^ annihilates each such vector. (C.f., 5 for more information.) As noted 
above, T-C™{n) is an invariant subspace for each n = 0, . . . , \ L/2\. 

One can define a subspace 'W^{n) to be the set of all vectors whose total 
8115(2) spin is L/2 — n. We would call this the "n-spin deviate" subspace because 
the total spin deviates from the maximum possible value of L/2 by n. (Total spin 
is a function of the Casimir operator, which generates the center of the algebra of 
SUg(2), which matches the usual notion for SU(2) total spin when q = 1.) Since 
the total spin operator commutes with S'^^ ^j, one can define subspaces H^'^{n, k) 
which are subspaces of 'W^in) with 5?^^ ^, eigenvalue equal to L/2 — k. These 
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subspaces are trivial unless n < k < L — n. Therefore 

L-2n 



k=0 

Also, W'^{n,n) = H'^^{n). For fc 1,...,L - 2n one has W'^{n,n + k) = 
(S'^ ^j)'^7i'™(n), and this is an isomorphic image. Since commutes with H^, 

each subspace is an invariant subspace for H^. 
It is natural to define 

E{L,n) = infspcc(i7'' f 7^"^(r^)) , 

which is the minimum energy of ranging over all vectors in the n-spin deviate 
subspace. We make explicit reference to the length of the chain [1, L] in this no- 
tation. On the other hand, commutes with and one can generate all of 
W'^(n) by acting on 7i'™(n) by some number of times. Therefore, it is clear 
that 

infspec {H^ \ = infspec {H^ \ H^'^in)) . 

This is the definition we will use henceforth. We now define = \ Ti}™{n). 
So, E{L,n) = infspec(iJ,'^). By Theorem 1.4 in [H], 

E{L, 0) < E{L, !)<...< E{L, [L/2\ ) , (13) 

for each finite L. By Proposition 4.1 in that same paper we know that all the 
inequalities are strict, at least as long as < g < 1. For q = Q one cannot define 
E{L, n) because there is no quantum group representation, but taking the limit as 
g — > 0+ gives E{L, 1) = • • • = E{L, [L/2\) — 1 which satisfy the inequalities, but 
not strictly. 

Moreover, by Proposition 7.1 in that paper, the sequence {E{L,n) : L > 
2n) is decreasing in L. Therefore, the limit lim l ~> oo E(L,n) necessarily exists. 
Obviously one has inequalities 

lim E{L,n + l) > lim E{L,n), 

L — >oo L — >oo 

which are derived from the fact that E{L, n + 1) > E{L, n) for every finite L. But 
note that one cannot automatically conclude that the inequality is strict in the 
limit. Whether or not this is so is a natural question. One might hope to resolve 
this question by finding an explicit formula for the limits. This is the first main 
result of the paper. 

THEOREM 2.1. For all n e N, andQ<q<l 

lim EiL^n)^^'-'^:f-'^^\. (14) 



6 



For two values of n a formula for Iuhl^oo E{L,n) was previously known. 
For n = 0, one obtains the ground state energy. It is well-known that the ground 
state energy is E{L,0) = for all finite L, in fact the Hamiltonian was con- 
structed to satisfy this condition. For n = 1, E{L,n) measures the spectral gap. 
The formula for this was calculated for all finite L in The value of the limit 
is limi^oo E{L, 1) = 1 — A~-^. Since A — ^{q + q^^), this is easily seen to agree 
with the result of the theorem. 

REMARK 2.2. A different method for bounding E{1) = limi^oo -E(L, 1) was 
given in |12j. That method is based on the martingale method, which proved 
useful in interacting particle systems Some inequalities of were made 
sharper in JTj. In particular, this led to an independent derivation of £'(1). 

REMARK 2.3. Let Eq{L,n) and Eq{n) = Ivchl^oo Eq{L,n) be the relevant q- 
dependent quantities for q £ [0, 1]. Using properties of the functions Eq{L,n) (in 
particular monotonicity in q, c.f.. Remark 15. 3() one finds 

lim lim Er{L,n) = lim Eq{L,n) < lim lim Er{L,n), 

r — >g+ L — >oo L — >oo r — L — ^cxd 

whenever the relevant limits exist. This can be used to recover the (obvious) fact 
that Eq{n) — 1 for q — Q and n > 1. Also it can be used to obtain the upper bound 
Eq{n) < for g = 1, which matches the obvious lower bound. For q = 1 direct 
spin- wave trial functions can also be used to verify limi_,oo Eq{L, n) = directly. 
(In fact, in 8 the finite-size scaling was calculated for Eq{L,n) when q = 1 and 
n = 1.) In our proof, we use the fact that < q < 1. Therefore, we will make this 
assumption, henceforth. 

2.2 Droplet Hamiltonians 

In |15j . two of the authors investigated low-energy vectors for three different 
droplet-type Hamiltonians based on the XXZ model. For (5 G M one may define 

Km ^ H^^M + ^ (1 - 5f - ^i) , (15) 

on W[i,L]- We call this the "droplet" Hamiltonian. Also on ^-[i^l] the spin chain 
with periodic boundary conditions (spin ring) is defined as 

Km = ^[1.^1 + ^^^^ ■ (16) 

We call this the "cyclic" Hamiltonian. Neither of these Hamiltonians has the full 
SUq(2) symmetry. But they both have the symmetry of S'^^ ^j. We define 

£^(L,n) := infspec(i/fi,ij fH[i.L](n)) and (17) 
E^y^{L, n) := infspec {HI^^^ \ H[i.L](n)) , (18) 
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for each L > n. 

Finally, let us recall that one of the ground-states of the infinite-volume 
Hamiltonian is the all-up-spin state. The GNS representation for this state will be 
constructed in Section|31 For now, let us write (TYz, wz, Hz) for the GNS represen- 
tation. There are subspaces Ti.z{n) which are invariant for the Hamiltonian, and 
such that Hz = 0^o'^^(")- These are the n-magnon subspaces. We define 

Ez{n) := infspec(iJz \ Hzin)) , (19) 

for each n. 

The second main result of the paper is the following. 
THEOREM 2.4. For S > 1, 

lim E\L,n) = lim E^^^L^n) = Ez{n) = lim E{L,n) . (20) 

L— i-oo L—^Qo L-H-oo 

3 Droplet energies in the infinite chain 

In this section we will give a precise definition for Ez{n), and calculate it for all 
n e N. 

3.1 Set-up 

We will start by constructing the GNS Hilbert space for the all-up-spin ground 
state of the infinite XXZ chain. Instead of following the usual GNS construc- 
tion, for this special case one can define the representation directly. We define 
the GNS Hilbert space, TCz as the direct sum of Hilbert spaces Hz{n) for n G 
N = {0,1,2,...}. Each of these subspaces are f^-spaces on countable sets. Let 
^0 = {0}- Define Hz{0) = ^^(A'o), which is a 1-dimensional space. For n G N+, 
let <Y„ = {a; G Z" : xi < ■ ■ ■ < Xn}, and define ?iz(n) = £^{Xn)- This defines 
Hz- Given n G N and x G Xn, define Sx & (-"^{Xn) so that 5x{y) — ^x,y These 
define the natural orthonormal basis. Note that the basis for Hz{0) is denoted ^0. 
Physically, the Xi are just the positions of the down spins. 

We should next define operators on Hz satisfying the same commutation 
relations as the spin- matrices 5* for i = 1,2,3 from the last section, except now 
for all X G Z, not just a finite set x G {1, . . . , L}. (These generate the C*-algebra 
on which the infinite XXZ Hamiltonian operates by the Heisenberg dynamics.) For 
each a; G Z, there is a representation of SU(2) on Hz, given as follows. For each 
n G N, if a; G Xn, then: 

• S^6x = m{x,x)6x, where m{x,x) equals +^ if a: G {xi, . . . , Xn} and — ^ 
otherwise; 

• S~6x equals if a: G {xi, . . . ,Xn}, and otherwise it equals (5^(2, .a,) where 
y{x, x) G Xn+i is 

yix,x) = {xi,...,Xk,x,Xk+i,...,x„), 
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for that k G {1, . . . , n} such that Xk < x < Xk+i (considering xq — —oo and 

Xn+l = +00); 

• S^Sx equals unless x € {ii, . . . , a;„}, and in that case it equals Sz(x.x) 
where z{x,x) G Xn-i is 

Z{x,x) = {xi,...,Xk^i,Xk+l,...,Xn), 

for that fc G {1, . . . , n} such that x = Xk- 

This is similar to the Fock space representation of the CCR algebra, except that 
there is a restriction to have at most one particle per site. Therefore, this is some- 
times called the hard-core Bose gas. (It is also related to a Fock space representa- 
tion of the CAR algebra using the Jordan- Wigner transformation.) 

The cyclic GNS vector is the vacuum vector, 6^ G Hz{0)- Then we may define 
the GNS Hamiltonian as 

Hz = ^ hx^x+i , (21) 
xei, 

where the interactions have the same formula as in Q , but relative to the present 
representation. As in the finite case, each n-magnon subspace is an invariant sub- 
space for the Hamiltonian. It will be convenient to adopt a notation for the re- 
striction to the n-magnon subspace 

Hzin) := Hz \ Hiin) . (22) 

We define the droplet energies 

Ez{n) = infspec Hz{n) , (23) 

for each n G N. The main purpose of this section is to prove the following result. 

PROPOSITION 3.1. For each n G N, 

(l-g^)(l-g") 

"""^"^^ (TT^XTT^- ^''^ 

The result is trivial for 71 = 0. Henceforth, we will consider n > 0. Let 
ei , . . . , e„ be the coordinate unit vectors in Z" . Then Hz {n) can be represented 
as a discrete integral operator by a kernel 

Hz{n)f{x) - J2 Knix,y)f{y), 
for all X ^ Xn- The kernel is given by 

n 

K^{x,y) = Y,Ki'Hx,y), (25) 

k=0 
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where: ^ ^ 

for fc = l,...,n — 1, 

and ^ ^ 

This kernel is symmetric, corresponding to the fact that Hj^{n) is self-adjoint. The 
kernel can also be used to define operators on ^p(<Y„) for p other than 2. One 
important preliminary stop is to observe that Hz,n is bounded. 

LEMMA 3.2. \\H^{n)\\ < n{l + A-^). 

Proof. Let Hz{n;p) be the operator on £P{Xn) with the kernel Kn, as above. So 
Hzin) — Hzin; 2). One knows that 

\\Hz{n;oo)\\ ^ max V \Kn{x,y)\, 

X^Xn 

and 

\\H^{n-l)\\ < max V \Knix,y)\. 

Since the kernel is symmetric, these two numbers - the maximum column sum 
and maximum row sum - are equal. Both are bounded by n(l + A~^). (They 
are actually equal to it.) This follows by considering the maximum number of 
off-diagonal entries in any row or column, which is 2n, as well as the max;imum 
diagonal entry, n. Both occur when cci+i > Xi + I for all i = 1, . . . , n — 1. By the 
Riesz convexity theorem, (c.f., |18| . Section 5.1), this gives the stated upper bound 
for ||i?z(n;p)|| for all p G [1, oo], in particular p~ 2. □ 



3.2 Direct integral, Fourier decomposition 

Let us define d = ei + ■ ■ ■ + En G Z". Note that if the coordinates of a; G Z" are 
ordered as a;i < ■ • • < a;„, then the same is true for x + d. Therefore, the function 
t{x) = X + d defines a bijection on = {a; G Z" : xi < ■ ■ ■ < x„}. It is trivial 
to check that Kn{x,y) = i4r„(r(a;), r(y)) for all a;,y G Xn- Therefore, defining 
T : TLz{n) Tizin) such that Tdx — '5r(x)i h follows that T and Hz{n) commute. 
The shift operator T has absolutely continuous spectrum. The analogue of block- 
diagonalizing Hz{n) according to the eigenspaces of T is to make a direct-integral 
decomposition of Hz{n) using the usual Fourier transform with respect to T. We 
describe this in some detail, next. 
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Let yn = {x E Xn '■ ^ X ■ d < n — 1} . There is a natural identification 
of this as the quotient space <-f„/T. Namely, for every x € there is a unique 
y G [y„ and fc G Z such that x = r'^iy). We define [a;] to be this y. Let l{Xn) be 
the set of all sequences on Xn- The operator T extends naturally to this vector 
space. We let £o{Xn) be the set of T-invariant sequences / G £{Xn). Moreover, we 
define io{Xn) to be the Hilbert space of functions in ioi'^n) such that the following 
norm is finite: 

yeyn 

Let = R/27rZ be the unit circle. Let L'^{S^,el{Xn)) be the Hilbert space 
consisting of the set of functions $ : §^ ^ ^q(A'„) such that for each x G Xn the 
function $(-)(a;) is measurable on The norm is 

||$f = £^|$Wf ^, 

where $(^) e ioi'^n) so ||$(^)1| is the norm in i^Xn). There is an analogous 
Banach space Lp(§^,£q(A'„)) for each p g [1, oo] such that 

\\m = l wmr^, 

where ||^'(0)|| is still the norm in £g(A'„). We continue to denote ||<i?||2 by just 
||$||. Then L°°(§^, £o(A'„)) is a dense subspace of L'^{S^,£Q{Xn)), as in the case of 
finite-dimensional vector- valued functions. 

We define a map, G, which is an analogue of the Fourier series, by 

g : L'^{S\£l{Xn))^f{Xn) : Q^x) = / e'^-<* $(^)(a;) ^ . 

Jo ^'^ 

As in the case of the usual Fourier series, we define the map on a dense sub- 
space first, and will eventually extend to the full Hilbert space using the isometry 
property. (Wc remind the reader to think of ("^{Xn) as fimctions on £'^{yn), where 
Af„ = X Z to see the analogy with the Fourier series.) One can determine that 
for these functions, 

1 n-l »27r in 

ll^^f = - E / E + 2'rfc/n)(y) g-^-'^f •<^/" . 

This map has a nontrivial null-space (except in the case n = 1) because of 
the sum over k. But it is if we restrict to the closed subspace L^S^ , iQ{Xn)) 
of L'^{S^,£l{Xn)), defined as all those $ such that 

$(6* + 2T:/n){x) = g-^'^''^"*/" ^{e){x) , (26) 
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then e is a partial isometry between Ll{S'^ , el{Xn)) Ci L°°(§\ £§(A'„)) (with the 
L^(§^, £q(A'„)) norm) and its range in L^{Xn)- One can then extend Q to all of 
Ll{S^ , £l{Xn)) , to obtain an isometry with £^(A'„). 

This map is surjective. For example, given any x E X, defining 

one can easily check that — 6^- Recall that £'^{Xn) is 'Hz{n). Let us define the 
map F : Hz{n) L^(§^, £§(A:'„)) as the inverse of Q. 

In a situation such as this, it is usual to call the Fourier-type decomposition 
a direct-integral decomposition, and to write 

n^{n) - / iliXr.)-. 

(This matches the notation of 16 , section 13.16.) Actually, we have a slightly 
more involved situation because of the constraint H26|) . We could reduce to the 
usual situation by restricting attention to e [0, 27r/n), which is a fundamental 
domain for such $. However, for notational purposes which arise shortly, we prefer 
to keep the present convention, and merely remember that H26|) must be satisfied. 

We can define a family of bounded operators Hz{n;9) on ^q(A'„), such that 
for every / G Hzin), 

TH^in)f{e) - Hz{n;9)Tf{0), 

for all 9 € S^. This can be done precisely because Hz{n) commutes with T. One 
usually then says that Hz{n) is a decomposable operator, and writes 

Jgi 27r 

We have the same caveat about remembering (|26ll as before. Particularly, we should 
pay attention to the fact that the proposed Hz{n; 9) should preserve this property. 

The operators H^{n; 9) are easiest to express as discrete integral operators, 
as was the case for Hz{n) itself. For each 6* G S"'^, we define 

n 

i^„,«(a;,y) = ^A'J](a),y), (27) 

fe=0 

where: 

for fc = l,...,n — 1, 
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and 



In addition to Hz{n;d) it is useful to consider the operator on the vector space 
£{Xn) defined through this same kernel. We denote this operator as JC{d). It is 
easy to check that the kernel does map the vector space to itself (i.e., there are no 
divergences, because all sums are finite). It is also easy to see that £o{Xn) is an 
invariant subspace. The restriction of IC{d) to £q(A'„) is Hz{n;6). 

From the decomposition, it is clear that for any / G Hzin), we have 

{f,Hz{n)f)^ {Tfie),H^{n,0)Tf{9)) —. 

Also, the map 9 Hi{n]9) is norm continuous. Therefore, by the Rayleigh-Ritz 
variational principle, 

infspec Hz{n) = min infspec Hz{n;9) . 

(At this point, the reader may be concerned that because of H26f) this might 
not be true. We leave it as an easy exercise to check that infspec Hz{n;9) = 
infspec Hz{n; 9 + 2'Kk/n) for all 9 and that moreover, the equation above is true.) 
We can see that the minimum is attained for 6* = (as well as possibly other 
values as per the last comment). This follows from the observation that for any 
/ S i^Xn), one has 

(I/I, i7z(n, 0)1/1) < {f,H^{n,9)f) , 

for all 9 € S^, where |/| e £l{Xr,) is the function |/|(x) = |/(a;)|. (This, in turn, 
follows because the kernel has nonpositive signs for the off-diagonal entries.) 
From this, one knows that Proposition 13 . II will follow if we show that 

3.3 The Bethe ansatz 

We will now develop the simplest possible application of the Bethe ansatz. The 
following is a well-known result, which we include for completeness. 

LEMMA 3.3. Let C \ {0}. Suppose that ^ e (C^)" satisfies, for k = 

l,...,n-l, 

e'^^k + e-^'^C+i = 2A . (28) 
Define the function, G i{^n) cls 

n 

/€(^) = n^r- (29) 

k=l 
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Then f is an eigenvector of K,{9) with eigenvalue equal to 

Eit ^) - E (l - ^ [e'' 6 + C'] ) • (30) 
fc=i ^ 

Proof. The key is to consider an operator on £{Z") which restricts to JC{d) on 
£{Xn)- Define the foUowing kernel on Z", 

" 1 

Let C{9) be the operator with this kernel. This is easily related to the Laplacian 
on Z". In particular, if one defines £ £{17^) by the same formula as in (|29|l 
(except on all of Z"), one has 

C{e)F^ = E{i,6)F^. 
For A: = 1, . . . , n — 1, define a kernel on Z" by 

Let M.k,k+i{9) be the operator with this kernel. Then we claim the following is 
true. First, for any x,y £ Xn, 

n-1 

Kn,eix,y) = Ln^0{x,y) -^M'^'g+^{x,y) . 

fe=i 

Second, if y G and a? G Z" \ then 

n-1 

The reader can check both properties easily. Because of these two properties the 
following is a fact. Suppose that F G ^(Z") is an eigenvector of C{9) with eigenvalue 
E{9), and suppose that A4k,k+i{9)F = for all fc = 1, . . . , n — 1. Then defining 
/ G £{Xn) to be the restriction of F, one knows that / is an eigenvector of IC{9) 
with the same eigenvalue E{9). 

The condition Mk,k+i{9)F^ — 0, is called the "meeting condition" in the 
context of the Betlie ansatz. It is 

2A F{y) = e''F{y + e,-) + e-''F{y - e^+i) , 

for every y G Z" such that yk+i — Uk + ^- For F = F^, this is equivalent to 

OA ^J^fc tVk + i _ ie/:Vk + i i^Vk + i , -ietVk tVk 

Dividing by Cfe+i^i this is precisely the relation in H28|) . □ 
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The condition in (|28|l is the same as the hnear fractional relation 

1 



Let us define the matrix 



e Kk+i 



A 







2A - e^^^fc 



1 



-1 q + q 

Then the linear fr'actional relation is expressed as 

A 







1 


1 







More generally, suppose Vk £ \ {0} is a vector such that vl/vl = e'^^fe. Then, 
defining Vk+i — Avk, we see that Vk+i £ \ {0} and ffc+i/wf+i = e**^fe+i. The 
eigenpairs of A are 



= q, v+ = 



-111 

,1/2 



and 



A- = q- 



-1/2 



Therefore, the most general solution to the linear recurrence relation Ufc+i = Av^ 
for fc = 1, . . . , n — 1 is 

Therefore, the most general solution to the linear fractional recurrence relation 
(EHl is 

yl/2 „fe-l/2 I _-l/2„-fe+l/2 

iBt ^ z ' q ' + z ' q 

^'^ ^1/2 ^fe+l/2 _|_ ^-1/2 g-fe-l/2 ■ 

(For this, we have taken a — z^l'^q^l'^ and /3 = z^^/^g^^/^, which is allowed since 
the two variables a and /3 only amount to one independent quantity in the ratio.) 
Let us define 

^1/2 gm-1/2 _|_ ^-1/2 ^-m+1/2 
^m{z) — ^m+1/2 _|_ ^-1/2 q-m-1/2 ' 

for all z. Then another way to write the most general solution of (|28() is ^fc = 
2m(fc) (z); where m(A:) = fc — (n + l)/2 and z G C. If we wish to have T-invariance, 
then we require ■ ■ • = 1. Since it is more convenient to work with e'^^/j, this 
is rewritten as (e'^^i) • ■ • (e*^^„) = e™^. One easily sees that 



:i-m{z)^-m+i{z)---^m{z) = 



^1/2 -A/-1/2 _|_ ^-1/2 „M+l/2 



;jl/2 qAf+1/2 ^ 2-1/2 ^-7\/-l/2 

So the condition for T-invariance is that z solves 



^1/2 ^-M-1/2 ^ ^-1/2 qM+1/2 
^1/2 qM+1/2 _|_ 2-1/2 q-M-1/2 
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where M = {n — l)/2. One solution is ^; = e'®, where 9 = Q{q, n, 6) is 



/ 1 + o" 

e = 2tan-i f tan(n6>/2) 



For each choice of N2, . . . , E N = {1,2,3,...}, there is a unique point 
y & yn with — Uh-i = N). toi k = 2, . . . ,N, and this labels all possible points 
in Therefore, one has 



n / n 



Ili=k 



leyn N2,...,NneK k=2 \j=k / k=2 \ llj=fel''. 



yey, 

But, for 2 < fc < n, we have 



j=k 



/2 _ -m+1/21 



cos(e/2)[g'"-i/2 + g-™+i/2] +isin(e/2)[q'"-i'' -g 
cos(e/2)[^^+V2 +g-M-i/2] +^sin(e/2)[5^+i/2 _ ^-M+1/2] ' 



for M ^ (n — l)/2 and m = {2k — n — l)/2. This clearly has norm less than 1. So 
ll/^ll, relative to iQ{Xn), is finite. 
The eigenvalue is 

^ / 1 \ 

m=-M ^ ^ 

where ^; = e*®. But one can easily veriiy that 
Therefore, by a telescoping sum 



1 



l + g2 



1 



1 



1 1 



1 + 1 + g'-"2: 



Putting this all together, we obtain the following. 

LEMMA 3.4. Given n G N+, we set M = {n - l)/2. For € (-7r/n, 7r/n), we 
define 



e = 2tan" 



1 + g" 
1-q" 



tan(n^/2) e (— tt, tt) , 
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and for m = —Ai, —M + 1, . . . , M we define 

^m-l/2gie/2 ^ ^-m+l/2g-ie/2 
"™ ~ (^m+l/2gje/2 _|_ ^— m-l/2g-i0/2 ' 

T/iera setting 

t.k — e ^fe-(„+l)/2 , 

for k = \, . . . ,n, we have that f^{x) ~ Y[k=i ^fe* defines a (normalizable) eigen- 
vector of Hi{n;d). Its energy eigenvalue is equal to 

F (0) = {l~q'){l~q'n 

' (1 + g2)(l + g»e*e)(l + g»e-^e)) ' 

In particular, when 9 — 0, this gives the formula from \24 



REMARK 3.5. Parts of this lemma are standard. For example, the linear fractional 
transformation was solved by Babbitt and Gutkin in However, they did not 
consider the direct integral decomposition. Instead they considered "generalized 
eigenvectors" . Also, no reference is made to the exact formula for the energy. A 
more explicit formula is 



From this we can derive the well-known dispersion relation for the isotropic model, 
lim^ii En{0) = -^{1 — cos9). However, to prove that this is the minimum energy 
for IIx{ri\9) is beyond our calculations because we would have to obtain the full 
diagonalization of Hz{n;9), which wc have not done. For q = 1, this was done in 
the important work by Babbitt and Thomas 2 . (The generalization to other q 
was outlined in T.) All that we need to calculate is En{9 = 0), which we handle 
by a different technique, next. 

REMARK 3.6. In |2I], Yang and Yang considered the ferromagnetic XXZ model 
to complement their famous and important work on the antiferromagnetic XXZ 
model 120]. They derived the linear fractional recurrence relation for 9 = with 
respect to the problem of calculating E'^^'^in). However, they did not solve the 
linear fractional recurrence relation, although they did set up a graphical method of 
solution which allowed them to determine the important fact that lim„i_^+oo 5„j = 
q and lim,„^_oo S,„ = l/q. 

This would have given them the result that lim„^oo E'^^'^{n) is finite. (In fact 
it is a = (1 — (j'^)/(l + q^).) But there is an unfortunate typographical error in 
their paper. They mistyped the formula for the energy in the equation just before 
equation (24) of [2] (compare to their definition in equation (11) of (Sill)- In our 
notation, their error is equivalent to saying En{^, 9 = 0) is equal to 

n 

fc=i 



17 



In other words, they left off an important factor (2A)~^. For this reason, they 
determined that the energy of an "edge spin" is asymptotically equal to — (2A — 1) 
(with our notation) instead of the correct value, which is 0. As a consequence the 
droplet nature of these states was not recognized at the time. 



3.4 Positive eigenvectors are ground states 

In Lemma rOl setting 6 = gives = and therefore ^fc > for all k. Hence is 
a strictly positive eigenvector of Hz{n; 0). Moreover, we claim that infspcc Hz{n\ 0) 
is the eigenvalue of This would be enough to prove Proposition 13.11 We also 
know that c\ — Hi{n] 0) is positivity preserving (when considered as an operator on 
^^(3^n)) and bounded. Therefore, the proof is completed by applying the following 
theorem. 

THEOREM 3.7. Let y he a countable set. Suppose that A is a positivity preserv- 
ing, hounded, self-adjoint operator on£^{y). If A has a strictly positive eigenvector 
f, then the eigenvalue of f equals the spectral radius. 

Proof. Without loss of generality, assume y = N+. Let us denote H = ^^(N+). 
For each N e N+, let Pjv be the orthogonal projection from Ti onto Hn — 
e'^i{l,...,N}). Let us define 

An := PnAPn- 
Obviously this is positivity preserving. Let us also define 

fN ■■= PNf. 

which is a strictly positive vector in Tijv- Let E be the eigenvalue of /. Then 

PNAf = EfN. (31) 
Let us define P^ = 1 — Pjv, and let us define another nonnegative vector in TijV: 

Jn := PNAP^f. 

Since P'^^f = f - f^, an obvious bound is WJnW < \\A\\ ■ \\f - /Ar||. Using lEU, we 
have _ 

ANfN + fN^EfN. (32) 

Let tjv — (/wi/jv) > and rjy = H/atIP > 0. We know that rjy is an 
increasing sequence with limit equal to H/p. Consider the operator An, which is 
self-adjoint ouHn, defined by 

ANg= ANg + rj/- Jn {fN, g) + fn {In, g) 
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Note that each of the three summands is positivity preserving, since /jv and /jv 
are nonnegative vectors. On the other hand, 

An In = An Jn + In + — /n ■ 

TN 

Therefore, by equation (|32|) . we have that 

, tN 



AnIn ^ [E 

TN 



TN ■ 



In other words, /n is an eigenvector of ^at. Since /n has strictly positive compo- 
nents in Hat, the Perron- Frobenius theorem guarantees that the spectral radius of 
An equals E + {In /tn)- From this we determine that 

lim maxspec(AAr) — E, 

N^oo 

because tn T II/P: while 

0<t^< ll/ivll- 11/^11 < Pll • 11/11- 11/ -/atII, 

and ||/-/iv|U as oo. 

We claim that p{A) < lmiN~>oo p{An) by the variational principle. This 
would imply that p{A) < E. We already know that E is in the spectrum of 
A, so that if p{A) < E we have p{A) = E. So we just need to prove p{A) < 
limAT^oo piAN)- 

To begin with, note that by the (Rayleigh-Ritz) variational principle for self- 
adjoint operators, we have 

p{A) = snp{{g,Ag) : \\g\\^l}. (33) 

Therefore, it suffices to observe that for any g £ Ti.^ we have that g = limjv^oo PNg 
and ^ 

Acj — lim An Pn g ■ 

N~*oo 

The fact that g = limjv^oo PNg is the usual density result (which we have already 
implicitly used). Since limjv_+oo ||5 — Pn gW^ = 0, one has 

lim WAg-ANPNgW^ = hm (.g - Pjv .g, (g - Pjv g)) 

< \\Af lim \\g-PN9\\'' = 0, 

and the perturbation 

[An - AN\g = r'^ [Jn [In, 9^ + In (/at, 5)) , 
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is bounded, in norm, by \\An- — An\\ < 2||/Ar||/||/Ar||. This converges to zero for 
reasons we have explained before. Therefore, for any g satisfying \\g\\ — 1, 

{g,Ag) = lim ( Pjv.g, 3 ) < lim p(Ajv) , 

as was claimed. □ 



4 Droplet energies in the droplet Hamiltonian and 
cyclic chain 

Our next goal is to compare the energies (recall their definitions in H15lll8|l ) 
E{L, n), E^y^iL, n) and E^{L, n), for 5 > 1, to Ez{n). The simplest case is E^{L, n) 
for S > 1. The desired result follows by two applications of the Rayleigh-Ritz vari- 
ational principle. 

LEMMA 4.1. For any 5 eR, 

Ei{n) > limsnp E^ {L,n) . 

L — >oo 

Proof. Given a < b, both in Z, define Xn{[a, b]) = {x <E : a < xi < ■ ■ ■ < Xn < 
b}. By definition, 

E^{n) - inf{(/, H^f) : / G Hzin) , ||/|| - 1} . 

Since the functions with finite support are dense in TCz{n), and since Hz is bounded 
on T-Cz{n), one can replace this by 

E^{n) = inf mf{(/,i7z/) : f E f {X.,{[-L, L])) , ||/|| = 1} . 

Moreover, by translation-invariancc one can shift a function on [—L, L] to a func- 
tion on [2, 2L + 2]. Therefore, 

E^in) = hm inf{(/,i7z/) : f e f {X^{[2, L - 1])) , ||/|| = l} . 
Suppose / e i'^{X[2x-i],n and ||/|| = 1. Define V £ by 

l<xi<-<a:„<L 

Then it will be apparent that {ip,hx,x+iip) = {f ,hx,x+if) for {x,x + 1} C [1,^]. 
Moreover, it is apparent that (■0,S'fV') = 1/2 and that = 1/2 because 

/ vanishes if there is any down-spin at sites 1 or L. Therefore, ^ip^H^^ i]^) ~ 

(V', i?[i,L]'/') and this is equal to {f,Hzf). So 
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In fact, by exactly the same argument, 

if,H^f) > El,{n), 
for any V > L. From this, one sees that 

mf{(/,i/z/) : /ef2(%L_i],„), 11/11 = !}> sn^E\L',n). 

L'>L 

But taking the hmit as L — > oo, we obtain the result. □ 
LEMMA 4.2. For any 6>1, one has 

Ej,{n) < liminf£;'^(L,n). 

n— >oo 

Proof. Suppose that tp £ 'H{L,n) is any normalized vector, so that HV'H = 1- We 
can write this vector as 

l<Xl<---<X„<i 

where F : X^i^L],n — > C is such that 

||i^f:= X \Fix)\' = l. 

Then we can extend -F to a function on Xn by defining 

' \0 if a; e A'„\A'[i,i],„. 

It is easy to see, from the construction of the single site representations of SU(2) 
on Hi, that (/, /ia;,x+i/) = {i^, hx,x+ii^) for {x,x + 1} C It is also easy to 

see that (/, hx,x+if) = for {x, a; + 1} C Z \ [1, L], because, in this case, / has no 
down-spin at either x or a; + 1. Furthermore, one can easily see that 

{fM,if) = {^,ht^) and {f,hL,L+if) = {i>,h+^) , 

where /i+ = TTy{hx,y Py ) is the partial trace of the XXZ interaction with P+ = 
51 + 5^ the projection onto the up-spin vector. But one easily calculates 

Therefore, one sees that the energy of / relative to is exactly equal to 
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with 5=1. Since the boundary field operator is positive semi-definite, one obtains 

< (V',i^[l,L]V') , 

as long as (5 > 1. From this, it follows that 

By minimizing over ip, and using the Rayleigh-Ritz variational principle, one ob- 
tains the result. □ 

Combining these two lemmas proves the following. 

PROPOSITION 4.3. For 6 > 1, limL^oo {L,n) exists and equals Ez{n). □ 

Let us now consider E'^^^{L^n). There is a proof completely analogous to 
Lemma [4. II for the following. 

LEMMA 4.4. Ez{n) > \im snp E'^y'' {L , n) . 

L — >oo 

To prove the analogue of Lemma 14.21 requires a different argument. Given 
any x € [1, L], we define the projections and P,J on 7i[i.L] by 

^l^-Sl and P- = 1-P+. 
Suppose that [a, b] C [1, L]. Then we define 

x—a 

Let us also define P[a,b] = P^^^ ^ ^^[a t] • following result follows from Corollary 
4.3in[TSl. 

LEMMA 4.5. Let he a self-adjoint operator on 'ii[i.L] o,nd M = \\Kl — Hl\\- 
Given E < oo, suppose G '^^U vector with WipW = 1 and {ijj,KLip) < E. 

Given any interval [a,a + £ ~ 1] C [1,L\ and any £' < £, there is a subinterval 
[b,b + £' -1] C [a,a + £- 1] such that 

\\PlbM+e'-i]ilA\ > 1 - 

where 

_ 2{E + M) 
' " j[£/£'\ ■ 

Moreover, when e < 1 , defining I = [b,b + £' ~ 1], 

{Pi^.KLPxi^) < O'^Kl^P) + (Me -t- 2(A-i + 2M)^e{l - e)) . 
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We can now prove the following. 
LEMMA 4.6. hminf £"^"{1,71) > Ez{n). 

L — 'oo 

Proof. The result is trivial for n = 0. Let us assume n > 0. 

Given L £ N+, let i/i^ be a vector in H[i^L]{n) such that HV'H = 1 and 
{'iljL,H'j^''^JL) = E''^''{L,n). We wiU now apply Lemma 1131 We take I ^ L no 
that a — \ and a + £ — 1 = L. We take £' = 2n. We have M = 1 because 
Hl - Hl^"" = hi^L and = 1. By Lemma lOI we can take E = 2Ez{n) as 

long as L is large enough. Hereafter, we assume that L is large enough. Then we 
obtain 

^ AE^jn) + 2 

We obviously have lim^^oo e^, = since — 0{1/L). On the other hand, given 
any interval X d [1, L] with |I| ~ 2n, we have 

Pii^L = 0, 

since -tpL only has n downspins, and Pj projects onto vectors which have downspins 
at all 2n sites in I. Therefore, we obtain 

Since we have translation invariance, we can translate the interval T as long as 
we translate ip^ by the same amount. Then, without loss of generality, we may 
assume that 

Let ip']^ ~ -^[i'ri]^[L-(i+i LJ^-E'- This can be written as 

for some e P'{^\n+\,L-n\,n)- It is easily seen that ||/^|| = HV'lII and that 
{f'L,Hzf'j^) = {ip'j^, H^^'^ip'j^). Therefore, using the second half of Lemma [4.51 and 
the fact that lim^^oo = 0, we conclude that 

liminfi;'=y^(L,n) > Ezin) , 

L—>-oo 

as desired. □ 

Thus we have proved the following. 

PROPOSITION 4.7. For each n G N, limL^oo E'y''{L,n) exists and equals 
Ez{n). □ 

REMARK 4.8. Yang and Yang actually calculated hm^— 00 ^^'^^'^(i, using the 
Bethe ansatz in I21j . However, their proof was more involved than ours, because 
they used the Bethe ansatz for finite L and took the asymptotic limit. Also, as 
mentioned before, because of a typographical error, they miscalculated E'^'^'^{n). 
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5 Droplet energies in the kink Hamiltonian 



highest weight vectors are vectors with total spin equal to the total S''^-component. 
We denoted the subspace of highest weight vectors on the chain [1,L] with total 
spin L/2 — n by 7i^^{n). We also defined the energy E{L,n) to be the lowest 
energy of the kink Hamiltonian H^, which is the restriction of H^^ ^ to T-C^'"{n). We 
are interested in the limit, limL^oo E{L,n). Throughout this section, we require 
< g < 1; for q = 1 we can use the SU(2) symmetry to show instantly that 
limL^oo E{L, n) = 0. Here is the main result. 

PROPOSITION 5.1. For each 7i e N and < q < 1, 

hm E{L,n) ^ E^{n) = , ^ ^ \ • (34) 

L^co (1 + + q") 

The essential part of this statement relates the ground state energy of the 
kink Hamiltonian for highest weight vectors with the ground state energy of the 
translation invariant Hamiltonion in the n-magnon sector. A detailed description 
of this relation is a bit lengthy but the first steps are well-known and standard. 
First we recall the construction of a basis for the heightest weight vectors for the 
finite chain. Then, we extend this basis to infinite volume and define the action 
of the infinite volume kink Hamiltonian, ^ . We then construct a continuous 
bijection, i?, between the highest weight vectors and the n-magnon vectors with 
the property that ff^n ~ Hi ^R- The main proposition is proved in Section 

ESI 

5.1 Basis for the highest weight vectors 

We use the notation [2]g = q + q~^ . Let us start by introducing a set of "valid 
brackets" from which we then define the "generalized Hulthen bracket vectors". 
They were first introduced by Lieb and Temperley in [^. Let n e N be fixed 
throughout this section. A valid bracket b = ([a;i,yi], . . . , [a;n,yn]) is a collection 
of n brackets [xi, yi\ with 1 < Xi < yt < L for all i and yi < . . .yn- In addition, a 
valid bracket satisfies the three conditions: 

1. (Exclusion) For each x G [1,-^], let (ib(a;) = #{i : Xi = x or yi = x}. Then 
db{x) < 1 for all x G [l,L]. 

2. (Non-crossing) If yi > Xj for some j > i then Xi < Xj < yj < yi. 

3. (Non-spanning) For any i if there is some x such that Xi < x < yi, then 
dbix) = 1. 

We define |b| = n. The set of all valid brackets on the chain [1,L] is denoted by 
V([l, L],n). Now, given a valid bracket b € V'([l, L],n), we define the "generalized 
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Hulthen bracket vector" ■^(b) e as 

i=l 

l'fl")[i,L] is, of course, the aU up-spin vector in 7i[i,L] = (C^)®-^. Lieb and Temperley 
proved that for any n G [0, [i/2j], the set V{[l,L],n) forms a basis spanning 
7Y'^'*^([1, L], n). Therefore, to calculate E{L,n), it suffices to find the minimum 
eigenvalue for the matrix of H^^ ^ in the basis of tp{b). In addition, the action 

of upon these basis vectors is simple and has a very appealing graphical 

representation which we will see after the following Lemma. 

LEMMA 5.2. For xg[1,L-1] and b e V([l, L],n), let 

Then (f> has the following values, depending on the case. 

• Ifdbix) = dbix + 1) = then (j) = 0. 

• If di,{x) ~ 1 but db{x + 1) = then there is a bracket in b, [xi,yi\, with 
Ui = X. Then </> = i^i^') where b' is defined relative to b by the replacement 
[xi,yi] -^[x,x + l]. 

• // db{x) = but db{x + 1) = 1 then there is a bracket in b, [xi,yi\, with 
Xi = x+\. Then (j) = tp{b') where b' is defined relative to b by the replacement 

[xi,yi] [x,x + l]. 

• If dbix) = db{x + \) = 1 then one of four possibilities occurs. 

— There are two brackets in b, [xi,yi] and [xj,yj], with yi = x and Xj = 
x+1. Then (p — i'ib'), where b' has the replacements {[xi, yi], [xj, yj]} — > 
{[x,x + l],[xi,yj]}. 

— There are two brackets in b, [xi,yi] and [xj.yj], with yi — x and yj = 
x+1. Then 4> = "^{b'), where b' has the replacements {[xi, yi], [xj, yj]} — *■ 

{[X,X + 1], [Xj,Xi\}. 

— There are two brackets in b, [xi,yi] and [xj,yj\, with Xi = x and xj = 
X+1. Then (f> = ip{b'), where b' has the replacements {[xi,yi], [xj,yj]} — > 
{[x,x + l], 

— The bracket [x.x + 1] is in b. Then (j) = —[2]qil){b). 

The proof is left to the reader who should notice that the vector {q~^/'^S~ — 
q^/'^Sy) |'ft')[i^i:,] is just a multiple of the singlet state between sites x and y. 

In the graphical representation proposed by Lieb and Temperley we consider 
the graph whose vertex set is the ordered chain [1,L\. The edges are the edges 
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specified by brackets. Namely, there is an edge between x and y if [x,y] G b. 
Because of rules (1-3) these are precisely the set of partial matchings which can 
all be made above the line passing through [1,L], and such that no two edges 
intersect, nor does any edge span an impaired vertex. We call the edges "arcs", 
and moreover since we specify that they should be above the line through [1,1/], 
we call them "upper arcs" . An example for L = 8 is 



{[3, 4] [2, 5], [7,8]} 




Then, if we define Ux,x+i = — [2]^ ^+1' '^^ have that Ux,x+i acts as follows. First 
of all, associate to each Ux,x+i a diagram of a pair of arcs; e.g, 





















1 i 


1 4 


1 i 


> « 


1 



Imagine concatenating this graph below the arc system for b. Then contract all 
loose ends down to vertices, and stretch the arcs to their normal shapes. Then one 
obtains the correct arc system corresponding to b'. The one exception is if b' = b 
because [a;, x + 1] S b. But in this case one obtains b with one additional circle, 
or "bubble". If one declares that the bubble takes a scalar value — [2]^ to remove, 
then one has the correct answer for (f) also in this case. 

^^^^ =-{<I + Q-')- 

At this point let us mention that the matrices Ux,x+i = — [2]g ^+1 satisfy 
the Temper ley-Lieb relations and that is an irreducible representation 

space of the Temperley-Lieb algebra. 
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REMARK 5.3. In this matrix representation, the only dependence on q conies from 
the "bubble" terms, which contribute a factor — [2]^ to the diagonal elements. Since 
this is clearly decreasing in q, it is trivial to see that the minimum energy of the 
kink Hamiltonian in any n-magnon subspace is decreasing in q. This is pertinent 
to Remark 12.31 

5.2 The kink Hamiltonian on Z in the Hulthen bracket basis 

The purpose of this section is to define the kink Hamiltonian on the Hulthen 
bracket basis vectors for brackets on Z instead of [1,L]. So let us first extend 
the notion of valid brackets and set V{[—L,L],n) to be the set of valid brackets 
b = ([zi, j/i], . . . , [a;„, y„]) on [— L,L] satisfying the same conditions (1~3) from 
above. Then we define the set of valid brackets on Z, V(Z, n) ~ Ulgn 

For any b £ V(Z,rt) we take some L so that [—L,L] contains all points in b 
and then unambiguously define the Hulthen bracket basis vectors ip{b) as in (|35|l 
and tensor with the all spin-up vector on Z \ [—L, L]. 

Now, for any x G [—L,L — 1] and b £ V{[—L,L],n) we define the matrix 
representation, 

b'eV([-L,L],n) 

Note that for any such lattice point x and bracket b we have the following fact: If 
L' > L, then by Lemma 

• If b' e V{[-L,L],n) then A|,(b',b) = Al{b',b). 

• If b' e V{[-L\L'],n)\V{[-L,L],n) then A|,(b',b) = 0. 
Thus we can define for any a; G Z the kernel by 

A^(b',b) = lim Al{b',b). 

L — *oo 

Moreover, by the same rules explained in Lemma 15.21 we have 

• If b e V([-i,L],n) and a; e Z \ [-L - then yl^(b',b) 
fa' e V(Z,7i); 

• If fa' e V{[-L,L],n) and x G Z \ [-i,L], then A^(b',fa) = 

fa G V(Z,7l). 

Therefore, we can define a valid, finite kernel, A : V(Z, n) x V(Z, n) 

A(b',fa) - E^"(&''&)' 

since all but a finite number of summands will be zero. Moreover, from this we 
conclude that if fa G V{[-L, L],n), then yl(fa', fa) = unless fa' G V([-i - 1, L],n). 



— for every 
= for every 
^ R by 
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Therefore, defining f(V(Z, n)) to be the set of all sequences on the countable set 
V(Z, n), we may define the linear transformations An and the infinite chain kink 
Hamiltonian ^ by 

fa'eV(Z,ri) 
b'eV(Z,n) 

Again, all but finitely many terms in the sum are zero. Since An and „ are 
defined by the same kernel they have the same spectrum. 

LEMMA 5.4. The linear operator An restricted to ^^(V(Z, n)) is bounded. 

Proof. For this we will use again the Riesz-convexity theorem. From Lemma |5. 21 
we observe the following facts: 

• If b = {[xi,yi], . . . , [xn, Vn]) 6 V(Z, n), then A'^{b', b) is zero unless |a; — < 
1 or |x — yi\ < 1 for some i G [1,?^]. In any of these cases, there is exactly 
one b' such that A^{b, b') is nonzero. 

• If b' = {[x[,y [],... jlx'j^^y'^]) G V(Z,?i), then A^(b',b) is zero unless x — x'^ 
and X + 1 = y'l for some i G [1, n]. In this case, the only way for A^(b', b) to 
be nonzero is if either: b — b'; or if b differs from b' by the replacement of the 
arc [x^,?/^] by another arc with one endpoint in {x^x + 1}; the replacement 
of the arc [a;-,?/-] and another arc, bracketing it, [x'^.y'j] by the arcs [x'j,x'j\ 
and [y[,y'j]. The total number of these possibilities is bounded by n. 

• Whenever A^{b, b') is nonzero, the absolute value is bounded by 1. 

Using these facts, one can conclude that An is bounded both in £^(V(Z,n)) and 
£°°(V(Z, n)) by bounding the maximum "row sum" and "column sum" ol A. There- 
fore, by the Riesz-convexity theorem, it happens that An is bounded on every 
£P(V(Z, n)) for 1 < p < oo. In particular, it works for p — 2. □ 

5.3 Bounds for the change-of-basis transformation 

We show now that the change from the Hulthcn bracket basis to the Ising basis in 
the 71-magnon sector TCz{n) is bijective and bounded. It is important that we do 
this on Z and to remember that q < 1. We use the notation from Section |2| 

LEMMA 5.5. Define a transformation R : £P{V{Z,n)) P{Xn) such that 
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Then R is invertible and bounded for p = 1 and p = oo. More precisely, 



Moreover, R intertwines between the operator An and Hx{n), namely, RAn = 
Hz{n)R. 

Proof. The intertwining property follows immediately since the kink boundary 
fields will telescope to in the i — > oo limit of the all-up-spin GNS representation. 

The colums of R represent the coefficients of ijj{b) in the expansion in the 
Ising basis. The map R is clearly injective. On the other hand, given an Ising 
vector Sx, we can write this as a norm-convergent telescoping sum of Hulthen 
bracket vectors. This is not possible on a finite chain nor when q = 1. 

Recall the following basic fact. If B : C"^ C"^ is a linear transformation, 

and 

ni 

{Bu)i = ^ BijUj , i = 1, . . . , n2 . 

Then 

\\B\\^ = max V \Bij\ . 

I<2<n2 ^ — ^ 

Furthermore, i < ||oo- [Although R itself is not a finite dimensional matrix, 
we could restrict the discussion of the changc-of-basis transformation to finite 
chains where we would show that our bounds are uniform in L.] 

If we expand 'tp{b) in the Ising basis Sx, then we see that the sum of the 
absolute values of the coefiicients is {q~^^^ + q^^^) , and we therefore have that 

\\R\\oo = {q-'/' + q'/T- 

The other bound is combinatorial and rests on the following claim. 

CLAIM: Let w = {wi, . . . ,Wn) be a collection of points in X^. There is an 
upper bound on the number of nonspanning systems of arcs 

ii = {[xi,yi],...,[xn,yn]), 

for which there exists a permutation n € Sn such that ■u'7r(fe) G {xk^Vk} for each 
k e[l,n]. The bound is (2n)!/n!. 

Given this estimate we instantly get 

,(2n)! 



\R\\i<\\R'U<q- 



Proof of CLAIM: Suppose that b = ([a;i,yi], . . . , [a;„,y„]) is a nonspanning 

system of arcs (i.e., satisfying conditions 1 and 3) but possibly with crossings, and 
such that there exists a n € Sn, such that Wk G {xT^^k), yT!-(k)} for all k € [l,n]. Let 
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Vk be the complementary point so that {wk^Vk} — {xTr{k),yir(k)} for each k. Let 
(vi, . . . , Vn) be the rearrangement of vi, . . . ,Vn in increasing order. Then we claim 
that the arc system b, whose arcs are (the ordered rearrangement of) {wk^Vk} is 
also nonspanning. The reason is that one can transpose the {vk,Vk+i) such that 
Vk > Vk+1, one-at-a-time to obtain the desired goal, and never span new sites. 
Particularly, the only new sites which could be spanned by transposing Vk > Vk+i 
are the sites in [vk,Vk+i]- But since Wk < Wk+i, these sites must have already 
been spanned by one (or both) of the two arcs whose endpoints are (the ordered 
rearrangements of ) {vk,Wk} and {vk+i,Wk+i}- 

Now, we can determine a set of points {u^ , . . . ,u^,Ui , . . . ,u^} such that 
{vi, . . . must be a subset of this one. This would obviously prove the claim, 
since there are at most (2n)!/(n!)^ such subsets, and b is uniquely determined 
by (wi, . . . ,Vn) which is obtained from (wi, . . . , w„) by permuting by one of the n\ 
permutations in Sn . 

Let {ui , . . . , u,^ } be the points such that is the first point to the left of 
Wn which is not among {wi, . . . , and for each k < n, the point u'^ is the first 
point to the left of Wk which is not among {wi, . . . , w„} or {w^^j^, ...,«,"}. We 
claim that if Vk < then it must be among {u^ , . . . , u^}. This is because, in this 
case, {vk, Wk} spans no sites other than {wi, . . . , w„} or {vi, . . . , Vn}- Since Vj < Vk 
for j < k, in fact it spans no sites other than {wi, . . . , Wn] and {vk, • . . , Vn\- A 
similar construction for {u^ , . . . allows the conclusion that if Vk > Wk then 
Vk is in the set {m^, . . . , u^}. 

□ 

REMARKS 5.6. 1. By the open mapping theorem, R^^ is also bounded. 
2. From the Riesz-convexity theorem we derive the bound 

\\R\\p<\\R\\\'''\\R\\]l' 

for any p £ [1, oo] with 1/p + l/q = 1. In particular, the map 

R:e{V{'L,n))^Hz{n) 
and its inverse are bounded. 

5.4 A Wielandt theorem 

This section extends a Wielandt-type theorem ^21 applicable to Banach spaces. 
We actually prove a stronger statement then needed. So let us consider a countable 
set X . Then, let /c : X x X — > M be a kernel with the following properties 

1. There exists a uniform fco < cxd so that < fc(x, y) < fco for all x, y G X; 

2. There is an integer N such that 

sup e X : k{x, y) ^ 0} <N, and sup #{y G X : k{x, y) ^ 0} <N . 
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Interpolating between £^ and £°°, we know from the Riesz-convexity theorem that 
this kernel defines a linear bounded operator, K : £'^{X) £'^{X). The first partial 
result concerns the spectral radius of a restriction of A'. So let F C X. Then we 
define the operator if to be the operator on £'^{Y) whose kernel is k\Y x Y. 

PROPOSITION 5.7 ( Generalized Wielandt theorem). We assume the same con- 
ditions on K as above. Let Y he a finite subset of X and let j > be a kernel on 
Y xY. Let J be the operator on £(Y) defined by this kernel. If j{x,y) < k{x,y) 
for all {x,y) GY xY, then p{J) < p{K). 

REMARK 5.8. The norm for £{Y) is immaterial since F is a finite set. 

Proof. By the standard Perron-Frobenius theorem for matrices, there is a vector 
ip G £{Y) with eigenvalue A = p{Y). 

By extending it to be zero on X \Y, we can also consider this as a vector 
in £'^{X). Moreover, by the our hypotheses, we have that Jip < Kij}. This implies 
that {K - X)^ > 0. Writing {K" - X")iP = K{K"-^ - X"-^)iP + X"-^{K - X)^ we 
conclude inductively that K^ip > X'^ip for all n G N. Since the kernel has positive 
entries, we get that \\K"'ip\\ > X"\\ip\\. Therefore, \\K"\\ > A", and 

p{K) = lim HAT" II > A = piJ). 

n — »-cxD 

□ 

REMARK 5.9. This proposition and proof follow 
5.5 Proof of the main proposition 

The last item of our business is to prove Proposition 15.11 about the ground state 
energies of the kink Hamiltonian. Since _ff„ is self-adjoint, it follows that 

infspec(iJz,r!) = inf {^,Hz,n'>P) . 

On the other hand, using Lemma 15.51 we can write 

infspec (iiz,n) = inf {(j),An4>) ■ 

0G^^(V(Z,n)) 
||i?0|| = l 

Moreover, with the natural identification of V([— L, L], n) C V(Z, n), we have den- 
sity 



£2(V(Z,n)) = cl U f{V{[-L,Lln)) 



\Len 



Therefore, 



infspec(-ffz n) = inf min {4>,A„ 

ieN 4,et.^(V([-L,L],n)) 
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Now, let e > 0. Since both Hz,n and An are bounded below, there does ex- 
ist an L G N and G P{V{[-L,L\,n)) such that \\R(j)\\ = 1 and {(j),An4>) < 
infspec(^„)-|-e. For such a vector (f> G £^(V([— I/, L],n)) we have a unique vector tp G 
H^^([-L,L],n) with the property that = 1 and = {-tp, Hf_^ ^^ip) . 

By shifting the interval [— L, L] to the right by i + 1 units, we conclude that 

E{2L + l,n) < e + infspec(iJz,«) • 

Since e was arbitrary, and using the monotonicity of E{L,n) in L, we have 

E{n) < infspec(ffz,n) • 

For the opposite inequality we use Proposition 15.71 Let X — V(Z, n). Let 
L G N, and Y = V{[—L, L],n). Then we consider the matrix B : Y xY M. given 
by the kernel 

L 

x=-L 

The operators J — n — B and K — n — An satisfy the conditions of Proposition 
15.71 Therefore, we conclude that 

p{n ~ B) < p{n - An) ■ 

But using the fact that the spectra of B and K are both real subsets (because the 
associated operators are similar to self-adjoint operators), we conclude that 

p{n — B) = n — infspec(-B) and p{n — K) = n — infspec(_R') . 

Therefore, 

E{L,n) > infspec(A„) , 

as desired. 

A Results for small q 

In this appendix, we collect some results for small q. The primary purpose of 
this is to verify the Bethe ansatz formulas for droplet eigenstates of the reduced 
Hamiltonian H{n,6) for other values of 6* G §^ than 9 = 0. Using the methods 
of Section 0] we can also treat the cyclic Hamiltonian H'^^'^ in the e** eigenspaces 
of the translation operator. In particular, the latter is interesting because this is 
the framework analyzed by Yang and Yang in '21'. Previously, this regime of the 
XXZ model (small q and cyclic boundary conditions) was rigorously analyzed by 
Kennncdy in using the methods developed in |7]. (The purpose of |S] was partly 
to give a pedagogic introduction to the methods of 7!, but it also gave new and 
interesting results for the XXZ model, some of which we describe below.) 
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Before going further, we would like to mention that in the paper proper, none 
of the arguments were perturbative. All applied to the entire region q G (0,1), 
which is the maximal interval where the results are valid. This is important to 
keep in mind when one considers the relatively simple arguments to follow. 

A.l Droplet energies in the infinite chain for small q 

Let us fix n G N+. Before stating the main result of this section, we recall the 
following. The Hamiltonians Hz{n,9) are periodic of period 27r/n in the sense 
that there is a unitary phase multiplication, as in H2t)|) . such that after conjugating 
by that Hi{n,6) and H{n,9 + 27r/n) are equal. In particular, this means that 
the spectrum is 27r/n periodic. Moreover, if there is an eigenvector of Hz{n,,9), 
then multiplying this eigenvector by the necessary phase produces the relevant 
eigenvector for H{n, 9 + 2-k jn). 

PROPOSITION A.l. There exists a constant — qo{n) > such that for < q < 
qo, the infspec of Hz(n,9) is an eigenvalue for all 9, and for 9 e (— tt/A^, tt/A^) the 
eigenvector is the one given in Lemma \'S.4\ The eigenvectors are norm continuous 
in 9, and are determined for all 9 using this and periodicity. Moreover, there is 
a constant 7(n, q) > such that there is a spectral gap above the ground state of 
Hi{n, 9) of size at least ^{n, q), uniformly in 9. 

Proof. Fixing n and 9, there is obviously a spectral gap above the bound state for 
the Ising model, g = 0. It is easy to see that the gap is 1 at g = 0. But the kernel 
Kn^e, when thought of as a function of q, varies in a way such that the associated 
operators are norm-continuous with respect to q, on £q(A'„). (As used before in 
the paper, this can be proved by obtaining row and column sum bounds, which 
pertain to £^ and £°°, and then using Riesz convexity.) Therefore, there is some go 
and some curve 7(71, q), positive for q < qg, such that IIz{n, 9) has a unique ground 
state and a spectral gap of size at least 7(n, q) for all 9 as long as < g < go- But 
the bound states found in Lemma var v continuously in g, therefore, they must 
be the actual eigenstates. □ 

REMARK A. 2. The argument of the proof is, to some extent, an analogue of 
Yang and Yang's argument from '51? but starting from the Ising model, not the 
XY model, and valid directly in the infinite volume limit. We would like to mention 
that more sophisticated and more powerful arguments of the Yang, Yang style were 
employed by Goldbaum in 0] to handle the more complex - but still Bethe ansatz 
solvable -Hubbard model. 

REMARK A. 3. Note that in the proposition above, one cannot choose go to be 
independent of n. The reason is that in our interpolation, the rowsums and 

columnsums do depend on n because of the number 0/ matrix entries. On the other 
hand, in ^Sj two of the authors proved a positive spectral gap for all g and n with 
g" small enough, which is uniform in this regime. Therefore, using that result and 
the present argument, one can obtain a single go which works for all n > 0. 
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Note that not only are the energies for the bound states continuous in 9 and 
q, they are easily seen to be analytic. This is simply because the kernel entries 
of the operator are analytic in 9 and q, and using the properties of the kernel 
(that there are a finite number of nonzero entries in each row and column) we 
deduce analyticity of (g-dependent) Fourier-reduced Hamiltonian Hq{n,9) in the 
weak-topology. Using the spectral gap this is sufficient to guarantee analyticity of 
the eigenvectors. Using analyticity in 9, we can obtain the following result. 

COROLLARY A. 4. The spectrum of Hz{n) in the range (0,ji{n,q)) is absolutely 
continuous. 

We will not give a detailed proof, but the reader is referred to Theorem 
XIII. 86 of USj. After conjugating by the spectral projection onto (0,7(ri, g)) the 
Hamiltonian satisfies the conditions of that theorem. 

REMARK A. 5. One probably expects that the entire spectrum of Hz{n) is ab- 
solutely continuous for all < g < 1. Using the results of 2 this is presumably 
provable at q = 1. But in general the translation-invariance suggests it is true. 

A. 2 The Hamiltonian for the cychc chain for small q 

Note that for the finite cyclic chain, just as for the infinite chain, there is a well- 
defined translation operator, commuting with the Hamiltonian. The following re- 
sult was proved by Kennedy in [Tj- 

PROPOSITION A. 6. There exists a qo such that for < q < qo the L lowest 
energy levels of H^^^ in the sector with n downspins (0 < n < L) can be indexed 

by the translation eigenvalues e**, for 9 — 27rfc/L and k G Z/LZ. For all 9, there 
is an analytic expression for the energy eigenvalue E'^y'^{L, n, 9) satisfying 



The coefficients dg — ds{n,q) are of order 0{q'^). 

REMARK A. 7. The arguments in 15^ prove that there is also a gap for large 
enough n and small enough g, and calculates the asymptotic form of the energy in 
the n ^ oo limit, with q fixed. As is easily seen from our present analysis, in that 
limit the energy converges to a = (1 — g^)/(l + g^). A simpler argument, but which 
is not robust to changes in L, can follow the proof of Proposition lA.il Namely, one 
can construct a kernel in each subspace of n downspins and translation eigenvalue 
0, and check that as a function of q the kernel is continuous, and moreover it is 
uniformly continuous for q G [0, 1]. If one considers the sequence of operators for 
all L (acting on different Hilbert spaces depending on L) one can even deduce that 
they are in some sense equicontinuous, because the stronger fact is true that the 
operators are Lipschitz with Lipschitz constants which are uniformly bounded in 
L G N+ and (? e [0,1]. 



DO 




s — — 



OC 
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REMARK A. 8. The important technique of Kennedy, which follows the previous 
work [7], is to obtain a perturbation expansion which can be performed for all L 
at once, therefore allowing comparison of different L. 

COROLLARY A. 9. The series expansion of Kennedy for the L — > cx) limit mat- 
ches the analytic expressions obtainable from Lemma \3.4\ 

Proof. One wants to show that 

lim E'y'{L,n,eL) = E^,„{e) . 

L — >oo 

One knows the existence of a spectral gap in the 9„ subspaces for small enough 
q uniform in L for H^{"]^] ■ In the last subsection, we established a similar spectral 
gap for Hz{n,0) on the infinite chain. Therefore, we can use exactly the same 
argument as in Section 0] to establish the same result for all G §^ that we 
established for 9 = 0, there: namely Proposition 14.71 The reader will find that 
translation invariance played no special role in that argument. □ 

Acknowledgement. B.N. acknowledges support and hospitality from the Erwin 
Schrodinger Institute for Mathematical Physics, Vienna, and the Centre de Phy- 
sique Theorique, Luminy, where part of this work was carried out. 

References 

[1] D. Babbitt and E. Gutkin. The Plancherel formula for the infinite XXZ 
Heisenberg spin chain. Lett. Math. Phys., 20(2):91-99, 1990. 

[2] D. Babbitt and L. Thomas. Ground state representation of the infinite one- 
dimensional Heisenberg ferromagnet. II: An explicit Plancherel formula. Com- 
mun. Math. Phys., 54(3):255-278, 1977. 

[3] L.D. Faddeev and L.A. Takhtajan. What is the spin of a spin wave? Phys. 
Lett, 85A:375-377, 1981. 

[4] P. Goldbaum. Existence of Solutions to the Bethe Ansatz Equations for the 
ID Hubbard Model: Finite Lattice and Thermodynamic Limit. Commun. 
Math. Phys., 258(2):317-337, 2005. 

[5] C. Kassel. Quantum Groups. Springer Verlag, 1995. 

[6] T. Kennedy. Expansions for Droplet States in the Ferromagnetic XXZ Heisen- 
berg Chain. Markov Processes and Rel. Fields, 11:223-236, 2005. 

[7] T. Kennedy and N. Datta. Instability of interfaces in the antifcrromagnetic 
XXZ chain at zero temperature. Commun. Math. Phys., 236(3):477~511, 2003. 



35 



[8] T. Koma and B. Nachtcrgaclc. The spectral gap of the ferromagnetic XXZ 
chain. Lett. Math. Phys., 40:1 16, 1997. 

[9] T. Koma and B. Nachtergaele. The complete set of ground states of the 
ferromagnetic XXZ chains. Adv. Theor. Math. Phys., 2:533-558, 1998. 

[10] S.L. Lu and H.T. Yau. Spectral gap and logarithmic Sobolev inequality for 
Kawasaki and Glauber dynamics. Commun. Math. Phys., 156:399-433, 1993. 

[11] T. Matsui. On ground states of the one- dimensional ferromagnetic XXZ 
model. Lett. Math. Phys., 37:397, 1996. 

[12] B. Nachtergaele. The spectral gap for some spin chains with discrete symme- 
try breaking. Commun. Math. Phys., 175:565-606, 1996. 

[13] B. Nachtergaele and L. Siegers. Construction of equilibrium states for one- 
dimensional classical lattice systems. // Nuovo Cimento, 100 B:757 778, 1987. 

[14] B. Nachtergaele, W. Spitzer, and S. Starr. Ferromagnetic ordering of energy 
levels. Joum. Stat. Phys., 116:719-738, 2004. 

[15] B. Nachtergaele and S. Starr. Droplet states in the XXZ Heisenberg model. 
Commun. Math. Phys., 218:569-607, 2001. 

[16] M. Reed and B. Simon. Methods of Modern Mathematical Physics, vol 4- 
Analysis of Operators. Academic Press, San Diego, CA, 1978. 

[17] W. Spitzer and S. Starr. Improved bounds on the spectral gap above frustra- 
tion free ground states of quantum spin chains. Lett. Math. Phys., 63:165-177, 
2003. 

[18] E. M. Stein and G. Weiss. Introduction to Fourier Analysis on Euclidean 
Spaces. Princeton University Press, Princeton, NJ, 1971. 

[19] H. N. V. Temperley and E. H. Lieb. Relations between the 'percolation' 
and 'colouring' problem and other graph-theoretical problems associated with 
regular planar lattices: some exact results for the 'percolation' problem. Proc. 
Roy. Soc, A322:252-280, 1971. 

[20] C.N. Yang and CP. Yang. One-dimensional Chain of Anisotropic Spin-Spin 
Interactions. I. Proof of Bethe's Hypothesis for Ground State in a Finite 
System. Phys. Rev, 150(l):321-327, 1966. 

[21] C.N. Yang and CP. Yang. One-dimensional Chain of Anisotropic Spin-Spin 
Interactions. III. Applications. Phys. Rev, 151(l):258-264, 1966. 



36 



